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Abstract 

In this paper, we consider the optimal disturbance rejection problem for possibly infinite 
dimensional linear time- varying (LTV) systems using a framework based on operator algebras 
of classes of bounded linear operators. This approach does not assume any state space repre- 
sentation and views LTV systems as causal operators. After reducing the problem to a shortest 
distance minimization in a space of bounded linear operators, duality theory is applied to show 
existence of optimal solutions, which satisfy a "time- varying" allpass or flatness condition. Un- 
der mild assumptions the optimal TV controller is shown to be essentially unique. Next, the 
concept of M-ideals of operators is used to show that the computation of time-varying (TV) 
controllers reduces to a search over compact TV Youla parameters. This involves the norm of 
a TV compact Hankel operator defined on the space of causal trace-class 2 operators and its 
maximal vectors. Moreover, an operator identity to compute the optimal TV Youla parameter 
is provided. These results are generalized to the mixed sensitivity problem for TV systems 
as well, where it is shown that the optimum is equal to the operator induced of a TV mixed 
Hankel- Toeplitz. The final outcome of the approach developed here is that it leads to two 
tractable finite dimensional convex optimizations producing estimates to the optimum within 
desired tolerances, and a method to compute optimal time-varying controllers. 
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Mathematical Preliminaries and Notation 



"B(E, F) denotes the space of bounded linear operators from a Banach space E to a Banach 
space F, endowed with the operator norm 

|| A|| := sup ||Ac||, A e B(£,F) 

xeE, \\x\\<i 

£ 2 denotes the usual Hilbert space of square summable sequences with the standard standard 
norm 



\x\ 



'■=^2\ x j\ 2 , x := (x ,x 1 ,x 2 ,- ■ -)e f 



3=0 

and inner product 



< x,y >=^2yiX h y = (2/0,2/1,2/2, ■ • -) G ^ 



Pk the usual truncation operator for some integer k, which sets all outputs after time k to 

P k x = (x Q , xi, x 2 , ■ ■ ■ , x k , 0, 0, • • •) 

An operator A e S(£ 2 ,£ 2 ) is said to be causal if it satisfies the operator equation: 

PkAP k = P k A, \/k positive integers 
and strictly causal if it satisfies 

Pk+\AP k = P k+1 A, , V7c positive integers (1) 

A partial isometry is an operator A on a Hilbert space which preserves the norm, i.e., ||Au|| = 
||w|| for all vectors u in the orthogonal complement of the null space of A [6J. 

A bounded linear functional / on a Banach space E is a real or complex valued linear map 
such that there exists c > with \ < f,x > \ < c\\x\\e, Vx G E, where < f,x > denotes the 
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image of x under /, and || • \\e is the norm on E. The dual pace of E, denoted E*, is the 
space of bounded linear functionals on E under the norm := sup a;g£ ; j » X » B | < /, x > \. 

Note that E* is itself a Banach space j6]. 

• A sequence {x n } in E is said to converge weakly to x G E if | < /, x n > — < f, x > | — > 0, 
V/ G E*. A sequence {f n } C E* is said to converge in the weak* topology to / G E* if 

| < f n , X > - < f, X > | > 0, Vx G -E. 

• A sequence of operators {T n } C -F) converges to an operator T G "B (E, F) in the operator 
topology if \\T n — T\\ — > 0. It is said to converge in the strong operator topology (SOT) 
if \\T n x — Tx||ir — y 0, Vx 6 £. K, denotes the space of compact operators in ¥>(E,F). An 
operator T G t B{E,F) is compact if there exists a sequence of finite rank (i.e., with finite 
dimensional ranges) operators {T n } C r B(E, F), such that, ||T — T n || — > as n — > oo [B|. 

The subscript " c " denotes the restriction of a subspace of operators to its intersection with causal 
(see [331 ED] f° r the definition) operators. "*" stands for the adjoint of an operator or the dual space 
of a Banach space depending on the context. denotes the orthogonal complement. The symbol 
"/\" denotes the closed linear span, and "/\" the intersection. denotes convergence in the 

weak* topology [6l [9] . Right and left hand sides are abbreviated as RHS and LHS, respectively. 

1 Introduction 

There have been numerous attempts in the literature to generalize ideas about robust control theory 
@H SB] to time-varying (TV) systems (for e.g. [121 H31 HH [33 HOI EHl EHl SOI EH and references 
therein). In [T^ Q2] and more recently [TT] the authors studied the optimal weighted sensitivity 
minimization problem, the two-block problem, and the model-matching problem for LTV systems 
using inner-outer factorization for positive operators. They obtained abstract solutions involving 
the computation of norms of certain operators, which is difficult since these are infinite dimen- 
sional problems. Moreover, no indication on how to compute optimal linear time-varying (LTV) 
controllers is provided. The operator theoretic methods in [T2l [T3~t [TT] are difficult to implement 
and do not provide algorithms which even compute approximately TV controllers. In [1] a state 
space extension of the Nehari Theorem to a time-varying system-theoretic setting is developed 
with parametrization of suboptimal solutions. In [10] the authors rely on state space techniques 
which lead to algorithms based on infinite dimensional operator inequalities. These methods lead 
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to suboptimal controllers but are difficult to solve and are restricted to finite dimensional systems. 
Moreover, they do not allow the degree of suboptimality to be estimated. An extension of these 
results to uncertain systems is reported in [H] relying on uniform stability concepts. 
In [8] both the sensitivity minimization problem in the presence of plant uncertainty, and robust 
stability for LTV systems in the £°° induced norm is considered. However, their methods could not 
be extended to the case of systems operating on finite energy signals. 

We believe that it is important to point to the fact that the lack of calculation of the degree of 
suboptimality in these methods, can result in arbitrary poor overall closed-loop performance. This 
fact is explored in [T8], where the authors, although in the linear time- invariant (LTI) case, studied a 
"two-arc" counter example that shows that in the limit suboptimal controllers can lead to arbitrary 
poor performance in the presence of plant uncertainty. It seems that this counter example has not 
yet received enough attention. 

Analysis of time- varying control strategies for optimal disturbance rejection for known time- invariant 
plants has been studied in [T6| |2] . A robust version of these problems was considered in [151 E] i n 
different induced norm topologies. All these references showed that for time-invariant nominal 
plants and weighting functions, time-varying control laws offer no advantage over time-invariant 
ones. 

In the book [35], the authors consider various interpolation problems in the stationary (LTI) and 
nonstationary (LTV) cases. In particular, Nevanlinna-Pick, Hermite-Fejer, Sarason, and Nudelman 
interpolations and Nehari extension problems are studied. In the TV case the unit disk is replaced 
by the set of diagonal matrices such as the weighted shift with spectral radius less than 1. The 
interpolation condition is replaced by an identity involving diagonal matrices. Two methods are 
used to solve these problems, a method based on the reduction to the time-invariant interpolation 
case but with operator valued functions, and another method based on a TV commutant lifting 
theorem developed by the authors of the book. The solutions for these problems are provided in 
terms of infinite dimensional operator identities. 

In this paper, we are interested in optimal disturbance rejection for (possibly infinite- dimensional, 
i.e., systems with an infinite number of states) LTV systems. These systems have been used as mod- 
els in computational linear algebra and in a variety of computational and communication networks 
as described in the excellent book [17J. This allows variable number of states which is predominant 
in networks which can switch on or off certain parts of the system [17], and infinite number of 
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states as in distributed parameter systems. We study questions such as existence and uniqueness of 
optimal LTV controllers under specific conditions. What is the corresponding notion for the Hankel 
operator which solves the optimal H°° problem? What does allpass mean for LTV systems? Is the 
optimal Youla parameter compact and if so under what conditions? Is it possible to compute op- 
timal LTV controllers within desired accuracy? Note that compactness is important since it allows 
approximations by finite rank operators (matrices), and therefore makes the problem amenable to 
numerical computation. 

Using inner-outer factorizations as defined in [21 E] with respect of the nest algebra of lower 
triangular (causal) bounded linear operators defined on £ 2 we show that the problem reduces to a 
distance minimization between a special operator and the nest algebra. The inner-outer factoriza- 
tion used here holds under weaker assumptions than [12], [13] , and in fact, as pointed in ([3] p. 180), 
is different from the factorization for positive operators used there. Duality structure and predual 
formulation of the problem showing existence of optimal LTV controllers is provided. The optimum 
is shown to satisfy a "TV" allpass condition quantified in the form of a partial isometry of an op- 
erator, therefore, generalizing the flatness or allpass condition [1H1 [2U ES]- The optimal controller 
and the corresponding dual operator are shown to be essentially unique under mild assumptions. 
With the use of M-ideals of operators, it is shown that the computation of time-varying (TV) con- 
trollers reduces to a search over compact TV Youla parameters. Furthermore, the optimum is shown 
to be equal to the norm of a compact time-varying Hankel operator defined on the space of causal 
Hilbert-Schmidt operators. The latter is a "natural" analogous to the Hankel operator used in the 
LTI case. An operator equation to compute the optimal TV Youla parameter is also derived. The 
results obtained here lead to a pair of dual finite dimensional convex optimizations which approach 
the real optimal disturbance rejection performance from both directions not only producing esti- 
mates within desired tolerances, but allowing the computation of optimal time-varying controllers. 
The numerical computation involve solving a semi-definite programming problem, and a search over 
lower triangular matrices. 

The results are generalized to the mixed sensitivity problem for TV systems as well, where it is 
shown that the optimum is equal to the operator induced of a TV mixed Hankel- Toeplitz operator 
generalizing analogous results known to hold in the LTI case [151 I3"8"l 121] . 

Our approach is purely input-output and does not use any state space realization, therefore the 
results derived here apply to infinite dimensional LTV systems, i.e., TV systems with an infinite 
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number of state variables. Although the theory is developed for causal stable system, it can be 
extended in a straightforward fashion to the unstable case using coprime factorization techniques 
for LTV systems discussed in [131 HI]- The framework developed can also be applied to other 
performance indexes, such as the optimal TV robust disturbance attenuation problem considered 
in [HI [19j [20] . For continuous LTV systems subject to time- varying unstructured uncertainty the 
problem was considered in [5], where it is shown that for causal LTV systems, it is equivalent to 
finding the smallest fixed point of a two-disc type optimization problem under TV feedback control 
laws. For the discrete version of the same problem, the duality structure of the problem and a solu- 
tion based on a new bilinear map was provided in terms of an infinite dimensional identity in [23J . 
A related problem involving the computation of the gap metric for LTV systems has been consid- 
ered in [22]. In particular, the computation of the gap metric for TV systems using a time- varying 
generalization of normalized coprime factorization as an iterative scheme involving the norm of an 
operator with a TV Hankel plus Toeplitz structure is provided. The duality structure is character- 
ized by computing the preannihilator and exploiting the particular structure of the problem. 
Part of the results presented here were announced in [311 [21] without proofs. 

The rest of the paper is organized as follows. Section [2] contains the problem formulation. The 
duality structure of the problem is worked out in section |3j In section HI the optimal solution is 
shown to satisfy a TV allpass property. Section [5] shows that under mild assumptions the TV opti- 
mal Q parameter is compact. Section [6] discusses the uniqueness of the optimal solution. In section 
El a solution based on a TV Hankel operator is derived together with an extremal identity for the 
optimum. Section [8] presents a numerical solution based on duality theory. A Generalization to the 
TV mixed sensitivity problem is carried out in section [HI Section [TO] contains concluding remarks. 

2 Problem formulation 

In this paper we consider the problem of optimizing performance for causal linear time varying 
systems. The standard block diagram for the optimal disturbance attenuation problem that is 
considered here is represented in Fig. [H where u represents the control inputs, y the measured 
outputs, z is the controlled output, w the exogenous perturbations. P denotes a causal stable linear 
time varying plant, and K denotes a time varying controller. The closed-loop transmission from w 
to z is denoted by T zw . Using the standard Youla parametrization of all stabilizing controllers the 
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Figure 1: Block Diagram for Disturbance Attenuation 
closed loop operator T zw can be written as [TT] . 

= Ti - T 2 QT 3 (2) 

where Ti, T 2 and T 3 are stable causal time-varying operators, that is, Ti, T 2 and T 3 e B c (£ 2 ,£ 2 ). 
In this paper we assume without loss of generality that P is stable, the Youla parameter Q := 
K(I + PK)~ l is then an operator belonging to B c (£ 2 ,£ 2 ), and is related univocally to the controller 
K [33]. Note that Q is allowed to be time-varying. If P is unstable it suffices to use the coprime 
factorization techniques in [391 [TT] which lead to similar results. The magnitude of the signals w 
and z is measured in the £ 2 -norm. Two problems are considered here optimal disturbance rejection 
which, and the mixed sensitivity problem for LTV systems which includes a robustness problem in 
the gap metric studied in [TTJ [36] . Note that for the latter problem P is assumed to be unstable 
and we have to use coprime factorizations. The performance index can be written in the following 
form 



/i := inf {||T ZW || : K being robustly stabilizing linear time — varying controller} 

inf ||Ti-T 2 QT 3 || (3) 

Q&h c (t 2 ,l 2 ) 

The performance index (J3]) will be transformed into a distance minimization between a certain 
operator and a subspace to be specified shortly. To this end, define a nest M as a family of closed 
subspaces of the Hilbert space I 2 containing {0} and i 2 which is closed under intersection and 
closed span. Let Q n := I — P n , for n = —1, 0, 1, • • ■ , where P_i := and P^ := /. Then Q n 
is a projection, and we can associate to it the following nest N := {Q n ^ 2 , n = — 1, 0, 1, • ■ ■ }. The 
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triangular or nest algebra T(N) is the set of all operators T such that TN C iV for every element 
N in J\f. That is 

T(A0 = {Ae-B{f,f):P n A{I-P n ) = 0,Vn} (4) 

Note that the Banach space H> c (£ 2 ,£ 2 ) is identical to the nest algebra T(J\T). Define 

AT = \j{N' e AT : N' < N}, N + = /\{N' e Af : N' > N} (5) 

where N' < N means N' C N, and N' > N means N' D N. The subspaces N N~ are called the 
atoms of M . Since in our case the atoms of M span £ 2 , then A/* is said to be atomic [3J. 
Following [3J we introduce inner-outer factorizations for the operators in T{M) as follows: 
An operator A in T{M) is called outer if the range projection, denoted P(Ra), Ra being the range of 
A and P the orthogonal projection onto Ra, commutes with A/* and AN is dense in the intersection 
N PI Ra for every N e A/". A partial isometry [/ is called inner in T{M) if [/*[/ commutes with A^ 
[H [3J [TT]- In our case, A e T(N) = , B c (£ 2 ,f) is outer if P commutes with each Q n and AQ n £ 2 is 
dense in Q n £ 2 (lA£ 2 . U € 23 c (£ 2 , £ 2 ) is inner if U is a partial isometry and U*U commutes with every 
Q n . Applying these notions to the time- varying operator T 2 G S c (£ 2 , £ 2 ), we get T 2 = T 2 iT 2o , where 
T 2i and T 2o are inner outer operators in ¥> c (£ 2 ,£ 2 ), respectively. Similarly, outer-inner factorization 
can be defined and the operator T3 = T^ T^, where T^i 6 S c (£ 2 , £ 2 ) is inner T 3o G r B c (£ 2 ,£ 2 ) is outer. 
The performance index [i in (J3j) can then be written as 

/i = inf || 7\ - T 2i T 2o QT 3o T 3i \\ (6) 

Following classical robust control theory (2H [251 HB] assume: (Al) that T 2o and T 3o are invertible 
both in £ C (£V 2 ). 

Assumption (Al) can be somewhat relaxed by assuming instead that the outer operators T 2o and 
T 3o are bounded below (see Lemma [I] p. 220). 

Assumption (Al) guarantees that the map Q — > T 2o 'B c (£ 2 , £ 2 )T 3co is bijective. In the time-invariant 
case this assumption means essentially that the outer factor of the plant P is invertible [25J . Under 
this assumption T 2i becomes an isometry and T 3i a co-isometry in which case T 2i T 2i = I and 
T 3 jTg i = /. By "absorbing" the operators T 2o and T 3co into the "free" operator Q, expression flS]) is 
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then equivalent to 



fji= inf II^TiT^-QH (7) 

Expression ([7]) is the distance from the operator T^TiT^ G !B(£ 2 ,£ 2 ) to the nest algebra 13 c (£ 2 ,£ 2 ). 
In the next section we study the distance minimization problem (J7|) in the context of M-ideals and 
the operator algebra setting discussed above. 



3 Duality 

Let X ba a Banach space and X* its dual space, i.e., the space of bounded linear functionals 
defined on X . For a subset J of X, the annihilator of J in X* is denoted J and is defined by [9], 
J 1 - := {$ G X* : $(/) = 0, / G J}, i.e., J is the set of bounded linear functionals on X which 
vanish on J. Similarly, if K is a subset of X* then the preannihilator of K in X is denoted K, 
and is defined by K := {x G X : = 0, $ G X}. The existence of a preannihilator implies 

that the following identity holds [9] 

min ||x — y\\ = sup | < x, k > | (8) 
y eK k&^K, \\k\\<i 

where < •, • > denotes the duality product. Let us apply these results to the problem given in (J7J) 
by putting 

X = !B(£ 2 , f), x = ItiTxT&i G S(£ 2 , f), J = S c (£ 2 ,£ 2 ) (9) 

Introduce the class of compact operators on £ 2 called the trace-class or Schatten 1-class, denoted 
C\, such that if {£;} is a basis of £ 2 , and T G Ci, then J2i < >< °°- Moreover, the sum 

< T^,^j > is independent of the choice of the basis [TU |3j. The operator T*T, T* being the 
adjoint operator of T, is positive (i.e., < Tx,x » 0, x G £ 2 and a; 7^ 0) and compact and has a 
unique positive square root, denoted (T*T)^. The trace-class norm is defined as [HE], 

llTUx := tr(T*T) L 2 = ^ < (T*T)^ i ^ i > (10) 
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where tr denotes the Trace. 

We identify 23 (£ 2 , £ 2 ) with the dual space of d , C*, under trace duality [14j, that is, every 
operator A in 23 (£ 2 , £ 2 ) induces a bounded linear functional on d as follows: $a e C\ is defined by 
<&a{T) = tr(AT), and we write T>(£ 2 , £ 2 ) ~ C* to express that "B(£ 2 , £ 2 ) is isometrically isomorphic 
to C\. 

Every trace-class operator T in turn induces a bounded linear functional on T>(£ 2 , £ 2 ), namely, 
$ T (A) = tr{AT) for all A in T>{£ 2 , £ 2 ). 

To compute the preannihilator of 23 c (£ 2 , £ 2 ) define the subspace S of d by 

5:= {Ted: (J - Qn)TQ n +i = 0, for all n} (11) 

In the following Lemma we show that S is the preannihilator of ¥> c (£ 2 , £ 2 ). 

Lemma 1 The preannihilator ofB c (£ 2 , £ 2 ) in C\, ±r B c (£ 2 , £ 2 ), is isometrically isomorphic to S. 

Proof. By Lemma 16.2 in J3] the preannihilator of T(W) is given by $t e 23 c (£ 2 , £ 2 ) if and only 
T belongs to the subspace 

{Ted : P(N~) X TP(N) = for all iV in e Af} 

where P(N) denotes the orthogonal projection on N, likewise for P(N~). P(N~) ± the complemen- 
tary projection of P{N~), that is, P(A^") ± = I - P(N~). 

In our case Af is atomic, and for any iV e Af there exists n such that iV = Q n+ i£ 2 , i.e., P(N) = Q n +i- 
The immediate predecessor of N, N~, is then given by Q n £ 2 , i-e., N~ = Q n £ 2 . The orthogonal com- 
plement of N~ is then £ 2 e QJ 2 - So P(iV-) ± = I - Q n . Therefore, 

S = {T e d : (/ - Qn)TQ n+1 =0,Vn} 

is isometrically isomorphic to the preannihilator of 23 c (£ 2 , £ 2 ), and the Lemma is proved. ■ 

The existence of a predual C\ and a preannihilator S implies the following Theorem which is a 
consequence of Theorem 2 in [9] (Chapter 5.8). 

Theorem 1 Under assumption (Al) there exists at least one optimal Q Q in T> c (£ 2 ,£ 2 ) achieving 
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optimal performance \i in fty, moreover the following identities hold 



fi= inf \\T*^T* ci - Q\\ = WT^T^ - Q a \\ = sup \tr{TT*fT^\ (12) 
Qez c {P,e?) TeS> ||t||i<i 

Theorem [1] not only shows the existence of an optimal LTV controller, but plays an important 
role in its computation by reducing the problem to primal and dual of finite dimensional convex 
optimizations. Under a certain condition in the next section it is shown that the supremum in (fl2l) 
is achieved. Theorem [1] also leads to a solution based on a Hankel type operator, which parallels 
the Hankel operator known in the H°° control theory. 

4 TV Allpass Property of the Optimum 

The dual space of S(£ 2 , £ 2 ), 3(£ 2 , t 2 )\ is given by the space, £(£ 2 , £ 2 )* ~ & ©i IC ± , where /C x 
is the annihilator of /C, the space of compact operators on £ 2 , and the symbol ©i means that if 
$ G C\ ©i K, 1 - then <3> has a unique decomposition as follows jT] 



$ = $ + $ r , ||$|| = ||$ || + ||$ r || (13) 

where $ G /C" 1 , and $ T is induced by the operator T e Ci, i.e, $ T (A) = tr(AT), A G S(£ 2 , £ 2 ). 
Banach space duality states that [9] 

inf \\x — v\\ — max |$(x)| (14) 
y&J $eJ x , ||*||<i 

In our case J = S c (£ 2 , £ 2 ). Since I>(£ 2 ,£ 2 )* contains Ci as a subspace, then !B C (£ 2 , £ 2 ) ± contains 
the preannihilator S, i.e., the following expression for the annihilator of J, J x , is deduced 

J 1 := 3 C (£ 2 , £ 2 )^ = 5©! (/Cfl S c (£ 2 , £ 2 ))^ (15) 

A result in [3] asserts that if a linear functional $ belongs to the annihilator J 1 - and $ decomposes as 
$ = $ D + $ T) where $ G (/C n S c (£ 2 , £ 2 )^ and $ T G S, then both $ G S c (£ 2 , i 2 ) L and $ T G 
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'B c (£ 2 , i 2 ) 1 - as well. We have then the following |3T] 




m II^T^-QII 



$ e (/C n 'B c ) ± , T e S 



max 



$o(T* T^) + tr{TT^T x T^ | (16) 



$o|| + II^Hl < 1 



If $ op4 = $ op £ i0 + ^T op t achieves the maximum in the RHS of (|T6|) . and Q the minimum on the 
LHS, then the alignment condition in the dual is immediately deduced from ( 1T61) 



Assume further that (A2): T^TiT^ is a compact operator. This is the case, for example, if T\ is 
compact, then ^^^(T^TxTg^) = and the maximum in ffTBl) is achieved on S, that is the supre- 
mum in (fT2l) becomes a maximum. Compact operators include, for example, systems with impulse 
responses that have measurable and square integrable kernels. 

It is instructive to note that in the LTI case assumption (A2) is the analogue of the assump- 
tion that T^TiT^ is the sum of two parts, one part continuous on the unit circle and the other in 
H°°, in which case the optimum is allpass [2"H [T8"] . We would like to find the allpass equivalent for 
the optimum in the linear time varying case. This may be formulated by noting that flatness or 
allpass condition in the LTI case means that the modulus of the optimum ((T^TxT^ — Q )(e id )\ is 
constant at almost all frequencies (equal to /i). In terms of operator theory, note that the optimum 
viewed as a multiplication operator acting on L 2 or H 2 , changes the norm of any function in L 2 or 
H 2 by multiplying it by a constant (=//). That is, the operator achieves its norm at every / e L 2 of 
unit L 2 -norm. This interpretation is generalized to the LTV case in the following Theorem, which 
part of it first appeared in j3T| without a proof. In fact, the optimal cost <^ is shown to be a partial 
isometry. 

Theorem 2 Under assumptions (Al) and (A2) there exists at least one optimal linear time varying 
Q e 3 C (£ 2 , £ 2 ) that satisfies the following 

i) the duality expression 



^(T'T^) + triT^T^T^)] = HT'T^ - Q \\ (||$ op , || + HT^y 



(17) 



max ||(/ -Q^T^Q 



(18) 
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holds, where T a is some operator in S and ||T ||i = 1- 

ii) and if \i > \i OQ := infgg^ HT^TiT^ — Q\, i.e. when the causality constraint is removed, the 
following operator 



1^21^1^301 Qo\ 



is a partial isometry. That is, the optimum is an isometry on the range space of the operator 
T Q in i). This is the time-varying counterpart of flatness. Note the condition fi > fi 00 is sharp, 
in the sense that if it does not hold there exist T^TiT^ and Q Q such that \^*^J^* ci _q° \\ ^ n °t 
a partial isometry. 

Proof, i) Identity (fTBl is implied by the previous argument that the supremum in (|T2l) is achieved 
by some T a in S with trace-class norm equal to 1. Combining this result with Corollary 16.8 in [3] 
(see also [Tj), which asserts that 

WT^T^ - Q a \\ = sup ||(J - Q^T^T+M (20) 



shows in fact that the supremum w.r.t. n is achieved proving that ( 1181) holds, 
ii) The operator T*T a is self-adjoint and compact, and therefore admits a spectral representation 
T*T D = Ylj ^jtftj ® 4>ji where {<pj} form an orthonormal basis for £ 2 consisting of eigenvectors of 
T*T D , and \j are its necessarily positive real eigenvalues [14J. It follows that 

E < T o T o<t>3, <t>i >= tr(fiT )$ = J2 X i = (21) 

3 j 

where we used the fact that T*T (pj = Xjcpj ^ 0. Let the polar decomposition of T Q be T Q = 
U(T*T )^, where U is an isometry on the set {4>j}, and (T*T )a the "square root" of T*T a . Now 
being the non-zero singular values of T Q , are also the eigenvectors of (T*T )2. It follows that 
for Xj ^ 0, we have 

1 -T (pj = -i=C/(T*T )5^. = J—U^/X^ = Ufa (22) 
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so {Utfij} is an orthonormal set which spans the range of T Q . Call ipj := U<f)j, then T a can be written 
as T = J2j ^j'Pj ® i>j and 



by definition of the trace 



< E ^ll^lhlK^i^M- Gorilla < E A ill^ T 3«-Q |lll^l|2 (24) 

< E A iH T i T ^-Qoll < \\T \\ 1 \\n i T l n ci -Q \\=^ (25) 



where the first inequality in ( 12411 follows from the Cauchy-Schwartz inequality, the second from a 
standard property of the induced norm, and the fact that {<pj} and {ipj} are orthonormal sets. 
Hence equality must hold throughout yielding 

im^T^ - q )^\\ 2 = nmti - Qowwu (26) 

for each ipj, that is, T^TiT^ — Q Q attains its norm on each ijjj, it must then attain its norm every- 
where on the span of {ipj}, so that [j^^przQ^ is an isometry on the range ofT a . 
Note that if condition /i > /i GO does not hold it is straightforward to find a counter example such that 
the optimal is not a partital isometry. Choose for example T^TiT^ to be a strict upper triangular 
which is not a partial isometry, then an optimal Q a is the zero matrix. The optimum is equal to 
llT^TiTg^H which does not correspond to a partial isometry. ■ 

Identity (Tl9|l represents the allpass condition in the time-varying case, since it corresponds to the 
allpass or flatness condition in the time-invariant case for the standard optimal H°° problem. In 
the next section, we show that the search over Q can be restricted to compact operators. This is 
achieved by introducing a Banach space notion known as M-ideals [3]. 

5 M-Ideals and Compact Youla TV Parameters 

In this section we show that the optimal TV Q is compact. This will allow approximations by 
finite rank operators to any desired accuracy, and consequently bring in finite dimensional convex 
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optimizations to carry out the computations. In order to achieve our objective we rely on the 
concept of M-ideals which is discussed next. 

Following [3] we say that a closed subspace M of a Banach space B is an M-ideal if there exists a 
linear projection II : B* — > M 1 - from the dual space of B, B* onto the annihilator of M, M 1 - in 
B*, such that for all b* G B*, we have 



||6*|| = ||n&*|| + ||6*-n6*|| (27) 

In this case M ± is called an L-summand of B*. The range N of (/ — II) is a complementary subspace 
of M ± , and B* = M ± (Bi N. A basic property of M-ideals is that they are proximinal, that is, for 
every b G B, there is an m in M such that inf mG jy/ ||& — m \\ — — m o\\- Under assumption (A2) we 
generalize Lemma 1.6. in [32] to causal LTV systems, i.e., the space H> c (£ 2 ,£ 2 ). Recall that Lemma 
1.6 states that if / is a function continuous on the unit circle, i.e., / G C, then 



inf \\f-g\\ 00 = M\\f-g\\ 00 (28) 

where A is the disk algebra, i.e., the space of analytic and continuous function on the unit disk 
A = H°° n C. That is it suffices to restrict the search to functions in A. To generalize (J2SJ) to causal 
LTV systems put B := JC, M := H c (£ 2 ,£ 2 ), and show that for b G K we have 

inf 116 — mil = inf 116 — mil (29) 

By Theorem 3.11 in [3j, S c (£ 2 , £ 2 ) n K is weak* dense in "B c (£ 2 ,£ 2 ). By Theorem 11.6 and Corollary 
11.7 in [3], S c (£ 2 , £ 2 )n/C, is an M-ideal in T> c (£ 2 , £ 2 ), and the quotient map q x : T> c (£ 2 , £ 2 )/(T> c (£ 2 , £ 2 )f] 
1C)i—> 'B C (£ 2 ,£ 2 )+JC/'B C (£ 2 ,£ 2 ) is isometric. Likewise, the quotient map q 2 : /C/(S C (£ 2 , £ 2 ) n/C)i — >• 
S c (£ 2 ,£ 2 ) + IC/'B C (£ 2 ,£ 2 ), is isometric. And the identity ( 129|) holds. In our case under assumption 
(A2) b = T*rr x T* ci G JC, and m = Q yields 

inf IIT*!^ - Q|| = inf H^TiI^ - Q|| = WT^T^ - Q \\ 

for some optimal Q Q G T> c (£ 2 ,£ 2 ) fl /C. That is under (A2) the optimal Q is compact, and thus it 
suffices to restrict the search in (13"U|) to causal and compact parameters Q. 

Now we turn our attention to studying questions regarding the uniqueness of the optimal solution. 
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6 On the Uniqueness of the Optimal TV Controller 



In this section under assumptions (Al), (A2), and fi > fi 00 , we prove that the optimal TV controller 
is essentially unique. Since there is an one-to-one onto correspondence between optimal controllers 
and Youla parameters, it suffices to show that the latter is essentially unique. Suppose by way of 
contradiction that there are two optimal compact Youla parameters Q\ and Q 2 such that 



/1 = wn^xTL - giii = wn^r^ - q 2 \ 



(30) 



and T^TiT^ — Q 1 and T^TiTg^ — Q\ are both isometries on the range of T , where T Q is the optimal 
dual operator in Theorem [2j In this case, any convex combination of Qi and Q 2 is also optimal 
since ||T*TiT 3 * { + AQi + (1 - X)Q 2 \\ > fi, VA G (0, 1), and 



\\\{T^ X T^ + (1 - A) (1*1^) + XQi + (1 - X)Q 2 \ 
< AH^T^ - Qx|| + (1 - A)||T 2 *T 1 T 3 * i - Q 2 \\ = p 



(31) 



implying that 



IT^T^ + XQi + (1 - A)Q 2 || = //, VA G [0, 1] 



In particular for Q 



* _ Q1+Q2 



, , let T := T 2 *TiT 3 * ci - Q*, and note that 
Q1-Q2 



r± 



By the parallelogram law we have for all x G 



(32) 



Q1-Q2 





2 






+ 

2 





Q1-Q2 



T n 3j 



2||r(T x)||2 + 2 



Q1-Q2 



(Tox) 



Since each of the operators T and T ± ^ 9 ® 2 is an isometry on the range of T a , we have 

= \\TT x\\ 2 = \\T x\\ 2 



r ± — - — T n x 



(33) 



(34) 



which implies that ||%^(T x)|| = 0, Vr G £ 2 , that is, we have necessarily Qi = Q 2} showing 
uniqueness on the range of T . ■ 
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Next, we show that the operator T Q in the dual maximization is unique. To see this, again suppose 
that there exists another operator T' Q such that 



\\T^T* ci - Q \\ = tr^T^T^ - Q )T )= tr^T^T^ - Q Q )Z) (35) 
By the proof of Theorem [2] there exist orthonormal sequences {(pi} and {4>j} such that 

t = y]\j<i)j®ij)j, \\t \\i = y]\j = l. (36) 

J 3 

Moreover since (T^T^ - Q )T = £\ ® {T^T*^ - Q )ipj, we have then tr^T^T^ - 
Qo)T )= J2j Xj < fa^T^T*^ - QMo >, and WT^T^ - Q \\ = // yields < fa, (T*!^ - 
Q )fa >= n, Vj. That is, (T*T 1 T 3 * . - Q )iffj = /ifa, Vj. The latter shows that (T^T*^ - 
Q )T = (T;T )K Similarly, (T*!^ - Q )T = (T>*T>)\ 

Now suppose that R := T Q — T' Q ^ 0, R G S 1 , i.e. i? is strictly causal, then ||T ||i = \\R + T^\\i = 
n = 1, and tr^T^T^ - Q ){T> + R))= tr ((T^T* ci - Q )T' ). Therefore, =\\T' + 

R\\ lf i = tr({T*fr x T* ci - Q )T' ), and then 

(T^T^-QjR = (T* l T 1 T.^-Q )(Z + R)-(T^T 1 T^-Q )Z (37) 
= {(T' + Ry{T' + R))^{TXY (38) 

showing that (T^TiT^ — Qo)-R is self-adjoint. But (T^TiT^ — Q Q )R is strictly-causal, that is, 
belongs to S, likewise for its adjoint {[T^TiT^ — Q )R) G S. This implies that 

((T*!^ - Q^RYdT^T^ - Q a )R) G S (39) 

Since ¥> c (£ 2 ,£ 2 ) is the annihilator of S it follows that 

tr^T^T^ - Q ) R) * ((T 2 * TiTg* j - Q )i?))= (40) 
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Thus, (T^TiT^ - Q )R = 0, and from ([38]) we have {{T' + R)*{T' + R)Y = {T' T' )* . Since 
T 2i TiT^ ci — Qo is an isometry on the ranges of T' Q and T Q = T' Q + R, we have for every i£f 2 , 

IKT^ + i?)^ = IKT^T^-Q^^ + ^xlh (41) 
= 11(25^ - Q )T' x\\ 2 = \\T' x\\ 2 (42) 

In the same manner, we get \\{T' a — R)x\\2 = \\T^x\\2- Next, apply the parallelogram law to get the 
equality \\{T' + R)x\\ 2 2 + \\{T' - R)x\\ 2 2 = 2||7>|| 2 , + 2||i2a;|||, from which we deduce that||i?x|| 2 = 
for all x G £ 2 , that is, we must have R = and then T = T' , showing uniqueness of T Q . ■ 

In the next section, we relate our problem to an LTV operator analogous to the Hankel opera- 
tor, which is known to solve the standard optimal H°° control problem [2"H HBJ- 

7 Triangular Projections and Hankel Forms 

Let C2 denote the class of compact operators on £ 2 called the Hilbert-Schmidt or Schatten 2-class 
[TU [3] under the norm, 

||A|| 2 := (tr(A*A)y, AeC 2 (43) 

Define the space A2 '■— C2 PI I> c (£ 2 ,£ 2 ), then A2 is the space of causal Hilbert-Schmidt operators. 
The orthogonal projection V of C2 onto A2 is the lower triangular truncation, and is analogous to 
the standard positive Riesz projection (for functions on the unit circle). Following [30] an operator 
X in B(£ 2 ,£ 2 ) determines a Hankel operator Hx on A2 if 

H X A = (I - V)XA, for A G A 2 (44) 

In the sequel we show that fi is equal to the norm of a particular LTV Hankel operator. To this 
end, we need first to characterize all atoms, denoted A n , of r B c (£ 2 1 £ 2 ) as A n := Q n +i — Q n , n = 
0, 1, 2, • • • . Write C + for the set of operators A G T> c (£ 2 , £ 2 ) for which A n AA n = 0, n = 0, 1, 2, • • • 
and let C 2 '■= C2 D C + . In [30] [29] it is shown that any operator A in 23 c (£ 2 , £ 2 ) D C\ admits a Riesz 
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factorization, that is, there exist operators A\ and A 2 in A 2 such that A factorizes as 

A = A 1 A 2l and ||A||i = ||Ai|| 2 ||A 2 || 2 (45) 

A Hankel form [■ , -] B associated to a bounded linear operator B G !B(£ 2 , £ 2 ) is defined by [301 [29] 

[A x , A 2 } B = tr(A 1 BA 2 ), A 1} A 2 G A 2 (46) 

Since any operator in the preannihilator S belongs also to ¥> c (£ 2 , £ 2 ) nCi, then any A G 5 factorizes 
as in (1451) . And if ||v4||i < 1, as on the RHS of ()12[) . A G S, the operators Ax and v4 2 both in A 2 
can be chosen such that ||Al|| 2 < 1, ||^4 2 || 2 < 1, and A n AiA n = 0, n — 0, 1, 2, • • •, that is, Ai G 
[29]. 

Now write V+ for the orthogonal projection with range A 2 C\C + . Introducing the notation (Bi, B 2 ) = 
tr{B\ Bi), and the Hankel form associated to B := T^TiT^, we have by a result in [29] . 

[Ai, A 2 ] B = tr{BA 2 A-,) = (A u (BA 2 y) = (V+ A h (BA 2 )*) = {A 1 ,V+{BA 2 )*) = (A±, H B A 2 ) (47) 

where H B is the Hankel operator (J — V)BV associated with B := T^TiT^. The Hankel operator 
H B belongs to the Banach space of bounded linear operators on C 2 . Furthermore, we have 



I^I^JI = sup \\H B A 2 \\ 2 = sup (A U H B A 2 ) (48) 

||A 2 ||2<l,A26Ca . A . 

|A 2 2 < 1,A 2 GC 2 



IA1II2 < l,Ax G C 



+ 



We have then the following Theorem which relates the optimal performance fi to the induced norm 
of the Hankel operator Ht^ jyt* .• The Theorem was announced in [31] without a proof. 

Theorem 3 Under assumptions (Al) and (A2) the following hold 

p = \\H T * TlT * || = || (/ - V^T.T^VW (49) 
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Proof. Since by the previous discussion any operator T G S can be factored as T = A\A 2 , where 
Ai G C£, A 2 G A 2 , \\Ai\\ 2 = \\A2W2 < 1, and ||T||i = ||Ai|| 2 || A 2 || 2 , the duality identity QE} yields 

fi = sup \tr(TT* i T 1 T* cl )\= sup HT^T*^^)! 

TeS, ||T||i<l 

P2II2 <i,A 2 e A 2 
\\Ai\\ 2 < 1,A 1 g C+ 

= \\Ht^t L X by 6HD 

By Theorem 2.1. [29] the Hankel operator is a compact operator if and only if B belongs to the space 
T> c (£ 2 , 1 2 ) + /C. It follows in our case that under assumption (A2) Ht*TiT* . is a compact operator on 
A 2 . A basic property of compact operators on Hilbert spaces is that they have maximizing vectors, 
that is, there exists an A G A 2 , \\A\\ 2 = 1 such that ||-^t 2 *TiT 3 * c J| = ll-^r 2 *TiT 3 * ci ^4||2, that is, a A 
achieves the norm of Ht*.t x t*^- We can then deduce from (fl2"j) an expression for the optimal TV 
Youla parameter Q as follows 

W H n^x nJ = \\H Ti . TlTic .A\\ 2 = \\H J * Tl7 * et _Q A\\ 2 = II (/ -V)(T* i T l T* ci A-Q A)\\ 2 
< WT^T^A - Q A\\ 2 < WT^T^ - Q \\ \\A\\ 2 < ||T*/riT* d - Q \\ = \\H T * iTlT *J 

All terms must be equal, and then 

|| (/ - V) (T^T* ci A - Q A) \\ 2 = \\T^T* ci A - Q A\\ 2 (50) 

Since T^T^A - Q Q A ={l-V) (T^T^A - Q Q A) +V (T^T^A - Q A) and V{T^{T^A - 
Q A)= 0. The TV optimal Q Q can then be computed from the following operator identity 

Q A = TZTxT^A - H T * iTlT ;A (51) 

The upshot of these methods is that they lead to the computation of fi within desired tolerances 
by solving two finite dimensional convex optimizations. 
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8 Numerical Computation of the Optimal Solution 

In this section a numerical solution based on duality theory is developed. If {e n : n — 0, 1, 2, • • • } is 
the standard orthonormal basis in £ 2 ,x then Q n £ 2 is the linear span of {e& : k — n + 1, n + 2, • • • }. 
The matrix representation of A e ¥> c (£ 2 ,£ 2 ) w.r.t. this basis is lower triangular. Note P n = 
I — Qn — > I as n — > oo in the strong operator topology (SOT). If we restrict the minimization 



in (H2D over Q G "B c (£ 2 , £ 2 ) to the span of {e n : n = 0, 1, 2, • • • , iV}, that is, Pj 



this yields 



a finite dimensional convex optimization problem in lower triangular matrices Qn of dimension N, 
that is, 



/ 



Q 



N 



Qn o 
Q21 Q22 








\ 



y <5ivi Qn2 Q 



N3 



Qnn J 



(52) 



and the optimization 



U.N '■= inf 

Qjve^^Pjv^.Pjv^ 2 ) 



T^iT^ - Q 



>N 



(53) 



where 



T 2i T l T 3ci -Q* 



( Tn - Q11 

T21 — Q21 
T31 — Q31 



TWl — Qjvi Tjv2 — QN2 Tn3 — Qn3 



T12 


T13 


Tin 


T 1(N+1) 


T22 — Q22 


T23 


T2N 


T2(N+1) 


T32 — Q32 


T33 — Q33 


T3N 


T3(N+1) 



l (N+l)l 



Ti 



(N+l)2 



T, 



(JV+1J3 



Tnn 



l N(N+l) 



T(JV+1)JV T( ]v+ i)( ]v+1 ) 



(54) 



where are fixed and correspond to the entries of T^TiT^, and Qij are variable. The optimization 
ix n overestimates /x, but fi^ I fx, and results in upper bounds and suboptimal TV parameters Qn,o 
and control laws, since Q is restricted to a proper subspace of r B c {£ 2 ) £ 2 ). By Arveson distance 
formula p] the minimization fl53|) is equal to 



Li N = max ||(/ - Q^iT^T^Q, 

Kn<N 



(55) 
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The degree of suboptimality can be computed explicitly as follows: Applying the same argument 
to the dual optimization on the RHS of ( fl2l) . by restricting S to the finite dimensional subspace 



S N := {T N E C x {P N t 2 , P N f) : (J - Q n )T N Q n+1 = 0, for all < n < N - 1} (56) 

In fact, with respect to the canonical basis {e; , i = 0,1, ■ ■ ■} for Pn^ 2 the subspace Sn is nothing 
but the space of N x N strictly lower triangular matrices. 
The dual optimization becomes 

fi' N := sup MT^T^M (57) 

where (T^TiT^J \ N is the restriction of the operator T^TiTg^ to PnC 2 - The supremum in f|5T|) is in 
fact a maximum, since (1571) is a maximization of a continuous linear functional over a compact set. 
The optimization ( 15?]) is a finite variable constrained convex optimization in the (strictly lowers 
entries of T/v), which yields lower bounds for /i since the dual optimization involves a supremum 
rather than an infimum, i.e., fi' N t /i, and suboptimal T/v i0 . 

The optimization problem fl57|) is in fact a semi-definite programming problem (SDP) since ||T/v||i < 
1 if and only if there exist matrices Y, Z 6 M Arx7V such that |12| 




try + trZ < 2 (58) 



The optimization f l57|) becomes then 

supremum tr{T N {T^T x Tl •) |^) (59) 
subject to fl58|) 

P n T N (I-P n+1 ) = 0, 71 = 0,1,-.. ,iV-l 

where P n is the N x N truncation matrix having as its first n columns the standards basis vectors 
{ej}, i — 1, 2, ■ • • , n and the remaining columns as iV x 1 zero vectors, that is, 

P n = (ei,e 2 , • • • ,e„,0, • • - ,0) 
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Under assumption (A2) T^TiT^ and the optimal Q a are compact, and by a compactness argument 
we have that Qn,o — > Qo in the operator topology, i.e., \\Q D — Qn,o\\ — > as N — > oo. Likewise 
Tn,o — > T Q as N — > in the trace class topology. Since Pjv — > I as N — > oo in the SOT. It is 
straightforward to show that these upper and lower bounds /ijy and fjf N converge to the optimum /x as 
N — > oo. These optimizations estimate // within known tolerance and compute the corresponding 
LTV operators Qn, which in turn allow the computation of LTV controllers K through the Youla 
parametrization. Solving such problems are then applications of finite variable convex programming 
techniques. 

Note the size of the truncations depends on the degree of accuracy desired. However, the size 
required for a particular accuracy can be exactly estimated by looking at the difference between 
the optimization ([55]) and the dual optimization (157|) . For periodic systems [ID], say of period q, 
it suffices to take the first q vectors {e n }^ =1 of the standard basis, i.e., N = q. In this case and 
for finite horizon problems the finite dimensional convex optimizations yield exactly the optimal 
corresponding TV Youla parameter Q Q and hence the optimal TV controller. 

9 The Mixed Sensitivity Problem for LTV Systems 

The mixed sensitivity problem for stable plants [HI HE] involves the sensitivity operator T\ := 

W \ ( W \ 

, the complementary sensitivity operator T 2 = P and T 3 := I which are all assumed 

o ; ^ v ) 

to belong to S c (£ 2 ,£ 2 x £ 2 ), and is given by the optimization 

= I l-l - \ PQ (60) 

where || • || stands for the operator norm in 23 (£ 2 , £ 2 x £ 2 ) given by 

||S||= sup (\\B lX \\l + \\B 2 x\\l)\ B= f Bl J (Gl) 
\\x\\ a <l, xeP x 7 \ B 2 

The optimization problem ( I6D1) can be expressed as a distance problem from the operator T\ to the 
subspace S = T 2 P S c (£ 2 ,£ 2 ) of S(£ 2 , I 2 x I 2 ). 

To ensure closedness of S, we assume that W*W + V*V > 0, i.e., W*W + V*V as an operator 
acting on £ 2 is a positive operator. Then there exists an outer spectral factorization Ai G S c (£ 2 , £ 2 ), 
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invertible in ¥> c (£ 2 ,£ 2 ) such that A}Ai = W*W + V*V pfl [TT]. Consequently, A X P as a bounded 
linear operator in ¥> c (£ 2 ,£ 2 ) has an inner-outer factorization UiG, where U\ is inner and G an outer 
operator defined on £ 2 [3]. 

Next we assume (A3) G is invertible, so U\ is unitary, and the operator G and its inverse G 
CB C (£ 2 ,£ 2 ). (A3) is satisfied when, for e.g., the outer factor of the plant is invertible. Let R = 
T 2 A^ 1 Ui, assumption (A3) implies that the operator R*R G H>(£ 2 ,£ 2 ) has a bounded inverse, this 
ensures closedness of S. According to Arveson (Corollary 2, [1]), the self-adjoint operator R*R has 
a spectral factorization of the form: R*R = A*A, where A, A -1 G ¥> c (£ 2 ,£ 2 ). Define R 2 = RA~ X , 
then R 2 R 2 = I, and S has the equivalent representation, S = R 2 H> C (£ 2 , £ 2 ). After "absorbing" A 
into the free parameter Q, the optimization problem ( IBTij) is then equivalent to: 

ti B = inf \\T X -R 2 Q\\ (62) 

Qe'Bcie 2 / 2 ) 

To solve the TV optimization ( IBIZj) it suffices to apply the duality results of section |3J The latter 
yields the predual space of S(£ 2 , £ 2 x £ 2 ) under trace duality, as the Banach space isometrically 
isomorphic to 



e 1 :={B=\ : Bt G CA (63) 




under the norm 



B\\ 2 U -^triB^ + B+B^ (64) 



The preannihilator ± S of 5 is characterized in the following Lemma. 
Lemma 2 The preannihilator ± S can be computed as 

L S = R 2 S® (I-R 2 R%)e 1 (65) 

where © denotes the direct sum of two subspaces. 

Proof. To show flS55 notice that tr($*T) = 0, VT in x 5 for $ G S(£ 2 , £ 2 x ^ 2 ), is equivalent to 
$*(! - = and = A* for some A G B c (£ 2 ,£ 2 ). This implies that $* = A*R*. By taking 

the adjoints we get $ = R 2 A G 5. 
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The following lemma is a consequence of (JH]). 

Lemma 3 Under assumption (A3) there exists at least one optimal TV operator Q Q G "B c (£ 2 ,£ 2 ) 
s.t. 

Ho = ||7i - R 2 Qo\\ = sup \trT*s\ (66) 

se x s,||s|[:i.i<i 

Note that Lemma yields [3] not only shows existence of an optimal TV controller K through the 
Youla parameter Q Q under assumption (A3), but also leads to two dual finite variable convex 
programming problems using the same argument as in section EJ which solutions yield the optimal 
Q within desired tolerance. A TV allpass property holds also for the TV mixed sensitivity problem 
under the assumption that T\ is a compact operator in 1> c (£ 2 ,£ 2 x £ 2 ) and fi Q > fi' OQ , where fj,' OD := 
infgg^ ||Ti — R2QW, i.e., when the causality condition of Q is removed. This generalizes the allpass 
property of section H] and is summarized in the following lemma. 

Lemma 4 If T\ is a compact operator and \x a > fi' oct then the optimal mixed sensitivity operator 
Tl ^ 2 ®° is a partial isometry from £ 2 onto £ 2 x £ 2 . The condition fi > /j,' 00 is sharp, in the sense 
that if it does not hold there exists T\ and R2 such that the optimum is not allpass. 

Proof. The proof is similar to the proof of Theorem 1 by working in the predual maximization. 
The details are omitted. If /i Q > ji' 00 does not hold, take for example V = 1 and W = ^=2, in which 
case Q = showing that the optimum is not allpass. 

Lemma H] is the TV version of the same notion known to hold in the LTI case as shown earlier 
in [27] using broadband theory and [47] . 



Next we show that the mixed sensitivity optimization is equal to the norm of a certain TV Hankel- 

R21 

Topelitz operator. In order to do this we use a standard trick in [llj. Since R2 = \ is an 



R 
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isometry, so is U :- 



R~2 
I — R 2 R 2 



. Thus, 



\i Q = inf \\U (7\ -R 2 Q) II = inf 



^ r* 21 w -q N 

y -R 22 R* 2l W j 



(I - R 21 R^ 2l )W 

Call n := I I . By analogy with the TV Hankel operator defined in section [3, we 

— R 22 R 21 W 

define the TV Toeplitz operator Tq*q associated to Q*Q as: 



T n * n A = Vn*nA, for AeA 2 



(67) 



The following Theorem solves the problem in terms of a TV mixed Hankel-Toeplitz operator. 
Theorem 4 Under assumption (A3) the following holds 



- \\ H R* iW H R * iW + Ti 



n*n\ 



(68) 



where Hr 21 + w is the Hankel operator associated to R 21 W , i.e., Hr* 2i wA = (I — V)R 21 WA, for A G 
A 2 . 

Proof. To prove ([6"8l define the following operator: 



T := IT 



W 




(69) 



where II is the orthogonal projection from A 2 x A 2 into the orthogonal complement of A 2 , A 2 x 
A 2 R 2 A 2 . The orthogonal projection II can be shown to be given explicitly by [23] 



n = I - R 2 VR 2 



(70) 



It follows from the Commutant Lifting Theorem for nest algebra (Theorem 20.22 p]) that 



fi = inf = \\T\ 



(71) 
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Since U is an isometry and thus preserves norms, operator pre-composition with U and using the 



explicit expression ( 1701) of II, straightforward computations show that \i 
therefore (I68D holds. ■ 



To 



and 



Theorem H] generalizes the solution of the mixed sensitivity problem in terms of a mixed Hankel- 
Toeplitz operator in the LTI case f27] 05j H7J to the TV case. This result also applies to solve the 
robustness problem of feedback systems in the gap metric [38] in the TV case as outlined in [Til EH] , 
since the latter was shown in [H] to be equivalent to a special version of the mixed sensitivity 
problem (I6TJ1) . 



10 Conclusion 

The optimal disturbance rejection problem for LTV systems involves solving a shortest distance 
minimization problem in the space of bounded linear operators. Dual and predual representations 
show existence of optimal TV controllers through the Youla parametrization. Under specific con- 
ditions, the optimal solution is compact and satisfies a "time- varying" allpass or flatness condition. 
The proposed approach leads "naturally" to a numerical solution based on finite variable convex 
programming, which involves solving a SDP problem for the predual and a search over lower trian- 
gular matrices for the dual problem. Moreover, a solution based on a TV compact Hankel operator 
is proposed. The latter leads to an operator identity for the optimal TV Youla parameter Q . A 
generalization to the mixed-sensitivity problem is carried out that draws certain analogies with the 
LTI case. 

Future work includes investigation of the numerical solutions by semi-define programming and find- 
ing efficient algorithms to take advantage of the special structure (lower triangular) of Q Q to solve 
numerically the dual problem. 
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